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Abstract 

We calculate the ground state energy of a massive scalar field in the 
background of a cosmic string of finite thickness (Gott-Hiscock metric). 
Using zeta functional regularization we discuss the renormalization and the 
relevant heat kernel coefficients in detail. The finite (non local) part of the 
ground state energy is calculated in 2+1 dimensions in the approximation 
of a small mass density of the string. By a numerical calculation it is shown 
to vanish as a function of the radius of the string and of the parameter £ 
of the nonconformal coupling. 



1 Introduction 

The Universe may have undergone a number of phase transitions since the big 
bang due to the spontaneous symmetry breaking in gauge theories. A number of 
topological objects may have been produced during the expansion of the Universe 
l] amongst which the cosmic strings seem to be of particular interest as seeds 
for the density fluctuations in the early Universe that are responsible for the 
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formation of galaxies (see for example ||). Also, the gravitational radiation 
produced by the formation of cosmic strings is part of cosmological scenarios. 

Several models of strings have been suggested. First of all Vilenkin has in- 
vestigated the case of an infinitely thin cosmic string 0. The energy momentum 
tensor of this string has a delta function like singular form. The space-time is 
locally flat except in the origin where it has a delta-shaped Riemann tensor |4]]. 
The vacuum expectation value of the stress-energy tensor for different kinds of 
fields in that background has been calculated both for zero temperature || and 
non-zero temperature || cases. 

The vanishing thickness of the string causes known problems. The vacuum 
expectation value of the stress - energy tensor has a non - integrabile singularity 
in the origin which can be seen already from dimensional considerations. As a 
consequence, for the renormalization of the ground state energy of quantum fields 
an additional counterterm is required. It is known as the topological Kac term 
0. This additional part may be recognized as due to boundary condition at the 
origin M. 

The mentioned problems can be avoided by considering a string with finite 
thickness. The simplest case is that of a constant matter density inside the 
string. It has been considered in Refs. || |§. The pressure p and energy density 
£ inside the string obey the condition p + £ = 0. The exterior of this string 
is a conical space-time and the interior is a constant curvature space ("cup" 
space). The metric is smoothly matched on the surface of the string but the 
scalar curvature has a laps on it. In fact, this space is a cone with a smoothed 
origin. There is no gravitational field outside the string in both above cases 
in opposite to the Newtonian logarithmic gravitational potential of a thread-like 
matter distribution. Note that this statement remains valid for an arbitrary radial 
matter distribution inside the string, provided that the translational invariance 
along the string is not broken. 

The purpose of this paper is to calculate the ground state energy of a massive 
scalar field in the background of a finite thickness cosmic string. 

In fact, we consider the (2 + 1) dimensional case. In zetafunctional regular- 
ization, the ground state energy of a scalar field $ is given by 

E = M 2 ° 1 -as- 1 -), (1) 

where 

c(*) = E(V)+™ 2 r (2) 

(«) 

is the zetafunction of the corresponding Laplace operator. The parameter M 
is arbitrary. It has the dimension of a mass. Usually it is denoted by which 
we reserve here for the linear mass density of the string. We assume the field 
$ to be put into a large volume with Dirichlet boundary conditions in order to 
render the eigenvalues discrete. It will be seen that the influence of this boundary 
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separates completely from the contributions of the pure background. The A( n ) 
are eigenvalues of the two dimensional Laplace operator 

(A - £R)y2 (n) (x) = A (n) (p(n){x). (3) 

where R is the curvature scalar. 

In the (3 + 1) dimensional case we would have to add the integration over the 
momentum of the translational invariant direction along the axis of the string 
and get 

^ 3+1) = ^ F(g ~ !\ m 2 *C(s - 1). (4) 
2 r(s-l/2) sv ; y 1 

The ultraviolet divergencies of the ground state energy are completely determined 
by the first few heat kernel coefficients. By means of 

oo , 

r dt t s 

o v ' 

and the asymptotic expansion for t — > 

K{t) = T^T E ^ (6) 
47rt n>0 

of the heat kernel K(t) corresponding to the operator in Eq. (||) the divergent part 
of the ground state energy can be expressed in terms of the first four coefficients 
(in (3+1) dimensions five coefficients would enter). We define 




+ B m wz Hi* ■ ( 7 ) 



In the following we calculate these coefficients for a string of zero thickness, 
reobtaining known results, and for a string of finite thickness in the approximation 
of a small mass density of the string. Using these coefficients, the renormalization 
of the ground state energy can be carried out in the standard way by adding the 
counterterms corresponding to (0) to a suitably defined classical energy. So we 
get the renormalized ground state energy 



77rcn 



lim (Eo - E^) . (8) 



In the (2 + 1) dimensional case we obtain the result Eq cu = in the given order 
of small mass density numerically. 

We note that the ground state energy defined in this way obeys the normal- 
ization condition 

^ cn ^0 for m^oo, (9) 
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which follows, from the circumstance that the heat kernel expansion is at once 
the asymptotic expansion for large mass. 

The organization of the papers as follows. In Sec.|2] we describe the Gott- 
Hiscock space-time of finite thickness cosmic string. In Sec.|3| we write down the 
general formulas for the thin string and calculate the corresponding heat kernel 
coefficients. In the next section we do the same for the finite thickness string. In 
Sec.|5| we calculate the ground state energy in the approximation of a small angle 
deficit. The result is discussed in Sec.|6] and an appendix contains some technical 
details of numerical calculations. 

We use units % = c = G = 1. 



2 The Space-Time 

The metric of Gott - Hiscock ||, |9| is a solution of Einstein's equations and it 
describes the space-time of an infinitely long straight cosmic string with non 
vanishing thickness. The energy density £ is constant inside the string and it is 
zero outside of it. The manifold can be covered by two maps. The first map, t G 
(— oo, +oo), p G [0, po],<p G [0, 2n], z G (— oo, +oo), covers the interior of the string 
and the second one, t G (— oo, +oo),r G [r ,+oo),<p G [0, 2ir],z G (— oo, +oo), 
covers the exterior. The coordinates (t, (p, z) are the same in both maps. r and po 
denote the radius of the string in external and internal coordinates, respectively. 
The string is situated along the z-axis. The metrics are C 1 - matched on the 
surface of the string, there is no surface stress energy (the extrinsic curvature 
tensors of the interior and exterior metrics are equal to each other f9j). The 
metric has the following form 



dsf n = dt 2 - dp 2 - p\ sm 2 ( — )dip 2 - dz 2 , (10) 

P* 



inside the string, and 
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ds 2 out = dt 2 - dr 2 - —dtp 2 - dz 2 , (11) 



v 2 



outside of it. Here p* = l/y8nS is 'energetic' radius of the string; £ is the energy 
density inside the string. The matching condition on the surface of the string 
links the exterior parameters (u, ro) and interior ones (p*,/?o) °f the string 

Po dcf , 1 r tane 

— = e = const , v = , — = . (12) 

p* cose po e 

From these relations we have the following consequences. The limit to the 
Minkowski space-time is achieved by letting the energy density inside the string 
tend to zero (p* — > oo) for fixed radius of the string p . Then the angle deficit will 
tend to zero too because of e — > 0. Thereby in this limit u = 1/ cose =1, ro = po 
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and both metrics turn into Minkowski space time. On the other hand side, in 
order to shrink the string (p — ► 0) at fixed exterior (e resp. the angle deficite are 
constant) we must turn the energy density £ to infinity proportional to e 2 /87rpg. 
Nevertheless the energy p per unit length of the string p which is the product of 
the energy density e 2 /Snpl and the cross section of the string is always constant 
and equals (1 — l/V)/4, the same value as for the infinitely thin cosmic string 
and it doesn't depend on the radius of the string. The two dimensional part 
{t = const, z = const) of the space-time (|T0| , |TTD is depicted in Fig.[l]. 



The manifold can be covered also by one map. One can continue the ex- 
terior radial coordinate r into the interior of the string by mapping r = r + 
(po/e) tan(ep/po — e) an d the space-time will be described by the metric 

d s 2 = dt 2 — — — dz 2 r G fO rnl 

[1 + e 2 (r - r ) 2 /p 2 } 2 v 2 1 + e 2 (r - r ) 2 /p 2 ' [U ' r ° J ' 

ds 2 t = dt 2 — dr 2 -dip 2 — dz 2 , r G [r , oo) . 

Here, the parameters u, po and r$ are connected by condition ([12]) . Nevertheless 
we shall use the metric in two maps because it is simpler for calculations. As far 
as the angle deficit is fixed, that is e = arccosl/V = po/p* is constant, one can 
exclude p* and rewrite the metric in the form which will be used in the following 

2 

ds 2 n = dt 2 -dp 2 - ^ sin 2 (— )d(p 2 -dz 2 , p G [0, p ] , (13) 
e 2 Po 



inside the string, and 



^2 

ds 2 t = dt 2 — dr 2 -dip 2 — dz 2 , r G [r , oo) (14) 

v 1 



outside of it. 



3 Infinitely thin cosmic string 

The metric for an infinitely thin cosmic string is given by (|14]) for r G [0, oo). In 
the (2+1) dimensional case, which we consider here, the coordinate z is absent. 
The eigenvalues A( n ) which enter the zeta function ([|) are determined by Eq. 
(HI) in the background of the spatial part of this metric. The curvature scalar 
R in this equation is proportional to the two dimensional delta function. So it 
is a potential with pointlike support and should be taken into account by a self 
adjoint extension of the operator corresponding to Eq. @. Instead, we drop 
these contributions here by considering the case £ = only. 
By means of the ansatz 

<P(n) = e m ^n(r) , (15) 
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constant curvature 
space - time 

— r = r , p = p 




r = R 



\ e = arccos(l/v) 



Figure 1: Plot of two-dimensional (t = const, z = const) part of finite thickness 
cosmic string space-time 



(n = 0, ±1, ±2, . . .) we arrive at the radial equation 

($_ _ n 2 v 2 - 1/4 A n _ Q 
\dr 2 r 2 J 

The solution regular at the origin of this equation is a Bessel function 



(16) 



n 



1 7r\r 



(17) 



We assume it to obey Dirichlet boundary conditions at r = R. Then the solutions 
A = \ n ,i (2 = 1,2,.. .) of the equation 



InXR 



J W (XR) = 



18) 



are the discrete eigenvalues A — >• A(„) = X n ^. Now the ground state energy and 
equally the zeta function (||D can be written in the form 



00 00 



(19) 



i=l n=— 00 

The sum over % can be rewritten as an integral (for details see |10| ) 





E' hin (s) = -M 2s 



cos [US 



2tt 



T,dnj m dk (k 2 -m 2 ) 1 / 2 -*— Ink-™ I nu (kR) , (20) 



where d n> o = 2 and do = 1 is the multiplicity of the angular momentum. Note 
the factor k~ nv in the argument of the logarithm, see [10 . 



In order to investigate the pole of i?o hm (s) we use the uniform asymptotic 
expansion of the modified Bessel function for n — > oo [ [TT|j 



t 



\nvz) 



nur/(z) 



2-nnv 



u k (t) 



k=l 



\nv) 



(21) 



with t = 1/Vl + z 2 , T](z) = + z 2 + \n(z/(l + y/1 + z 2 )) and z = kR/nv. 

The pole term of the zeta function ((s) (0) in s = —1/2 resp. of the energy 
Eq(s) in s = will be delivered by the first few terms (up to k — 2) of this 
expansion when inserting them into E (^0|). We note that this expansion is 
at once an asymptotic expansion for large masses and that the individual heat 
kernel coefficients enter expansions like (|7D multiplied by the corresponding power 
of m. Therefore we will in the following keep track of these powers and drop all 
contributions which for m — > oo are of order 0(— ). 

After inserting the expansion (^l|) into (l20f) , the integration over k can be 
carried out using 



dx{x 2 -l) 1 l 2 ~ s x{l + x 2 /a 2 )- p l 2 



r(| 



s)T(s 



2r(§; 



P~3 \ 



P-3 
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and we obtain the following expression for E(s) 



E(s) 
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r( s )z(o, s )-^^±^z(2, s + i) 
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884, 



25r(s + l)Z(0,s + 1) 



2)Z(2, s + 2) + — r(s + 3)Z(4, s + 3) 



r(s + 4)Z(6,s + 4) 



+ . . . 



(22) 



Here 2 F\ = 2 Fi 



i S - I" -■ 

2 ' > 2' 

/ ni/ \ 2\ 



is the hypergeometric function; Z(p, q) = 

, , , ~ q and 6 = mR. Next we have to perform the analytical 
continuation s — > 0. 

First of all let's consider the part containing the hypergeometric function 



rfs-r 
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+ 2E 



n=l 



r(^-i) „ ^ 1 

-2-^1 — r r ( s - t: 



(23) 



For the calculation of the series we use the Mellin 
of the hypergeometric function 



1) 



1 



1 1 

■ — , s — 1: -; 

2' ' 2' 



Barnes type representation 

tJ J 



1 fT(s-l + t) 
t - 1/2 



2vri 



r(-t) f 



2/ 



where the contour is such that the poles of T(s — 1 + t)/(t — 1/2) lie to the left of 
it and the poles of r(— t) to the right [fTH| . Before interchanging sum and integral 
one has to shift the contour 7 to the left crossing the pole at t = 1/2 up to 
t = —1/2 and then to close it to the left. This is because of the convergence of 



the series 



00 

n=l 
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it is necessary to have t < —1/2. The residue at the point 

) . Taking 



t = 1/2 cancels the second manifestly divergent term in the sum 
the limit s — > we obtain the following finite expression for (|23| ) 
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no) = -7= t c R (-2) + -s 



E 



1 / 2tt6 
-2) + — 2 In 

'7T \ ^ 



+ 



(-1)" Cg(2n + 2) 
0F^ O (n + l)(n + 2) n + 3/2 



2+2rt 



Here (n(a) is Riemann zeta function. 

The series is absolutely convergent for bjv < 1, but we need it in the domain 
of large R, that is for b ^> 1. For this reason we perform the analytic continuation 
into the domain we need. With t = ibj v in the formula 

00 j2n 



E — C*(2n) 



In 



irt 



n=l 



?? 



sin Tit 



\t\ < 1 



(24) 



we can express the series in Y(0) in terms of the series ([3). As a result, 3^(0) is 
expressed in terms of functions which are analytical in the whole plane of b and 
it can be divided into a polynomial and an exponentially small part 



Y(0) 
with 
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3b 



In (l - e~^) - 2v^F [4Qx(b/v) - 2Q 2 (b/u) 



Qa(x) = — 

x a Jo e 



dtt a 



2-rrt 



- 1 



1 T(a + 1) 1 

^72^ CMa + 1) - - ^ 



dtt a 

3 2yrt _ 



(25) 



For large 6 this expression is exponentially small. Therefor we arrive at 

b 



Y(0) 



3b 



0(e 



-b\ 



s 



The same result may be obtained in another way. One can use an analytical 
continuation of the hypergeometric function JET], namely 



+ 



2' ' 2' V b J J b T{s - 1, 

r(i/2)r(i/2-s) / (nu\*y- a ( 1 1_ 

r(-i/2)r(3/2 -s){ 1 + \TJ ) 2Fl | v 1 ' s ~ 1;s + 2 ; 77^) 



The first term in the rhs. cancels the second, divergent term in the sum 
Next, one can use power series expansion for the hypergeometric function because 
its argument 1/(1 + (nu/b) 2 ) is always smaller then unity The result will be the 
same as that obtained above by a longer calculation. 

The series Z(p, q) in (p2|) can be expressed in terms of the Epstein - Hurwitz 
zeta function 11121 



( E H(r,c) = Y,(n 2 + 



c 2 )' r 



n=l 



A known, quickly convergent expansion for large values of the parameter c is 



2 2r(r) r(r) ^ 



-2r 



0FT(r- 1/2) , r+1 , 27TV-+V2 « 



For small c it holds 



00 F(' r7 _|_ r ^ 

Ur.c) = £(-i) n 1 Vfr c2r ^ 2r + 2n ) • 

n=0 1 V) n - 

For integer r > 1 the zeta function can be expressed in terms of elementary 
functions. So we get the relevant part of the ground state energy in the form 

nthin / \ ( M \ 2S 1 ( 27ri? 2 , TT 3 / 2 Rm 2 r(s - 1) 77/ 1, 



m/ 87T y 3z/ z/T(s — |) 3 

7r 3 / 2 r(. 

32RvT(s 



+ ,o^ (a) iJ - ^(TT7)- ™ 



2 ' 
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dropping contributions proportional to exp(— Rm). By comparing this formula 
with Eq 1v fl7|) we can read off the heat kernel coefficients 

nR 2 7l 3 / 2 R 7T / 1\ 7T 3 / 2 , . 

Now, taking into account the general structure 

B r = [ c r dS + / a r dV (28) 

JdV JV 



of the coefficients we represent B\ as 



2tt 7T / 1 , 

Bx = — + - k-- . (29) 



3u 3 V z/ 

In this representation, the first term in B\ and all other coefficients are seen to 
be the result of the boundary at r — R. In fact, they are known fI5| . The second 
term in B\ is independent on that boundary. It is known as the topological Kac 
term M and is a result of the conical singularity 

Using these coefficients, by means of Eq.(|7]), we can define E$ lv and the renor- 
malized ground state energy E™ 11 (§). Now we observe, that all contributions in 
( |26|) except for the topological Kac term are due the the boundary at r — R. 
Leaving them aside, only the Kac term remains. It must be included into Eq 1v 
and we get 

^ cn = (30) 

for the genuine contributions of the string, or equally in the limit R — > oo. Let 
us remark, that this result holds also in (3+1) and higher dimensions as can be 
inferred from dimensional reasons. The parameters entering the problem are the 
mass density of the string which enters together with the graviational constant 
to form a dimensionless combination, expressed by the angle deficite, for example 
and the mass of the quantum field. In the case of a string with zero radius there is 
no further dimensional parameter on which the renormalized ground state energy 
might depend. As the ground state energy has the dimension of a inverse length 
it might be proportional to the mass but such terms had been subtracted in the 
renormalization process. Note that this discussion does not apply to the case 
of a string with finite radius because this radius is the additional dimensional 
parameter allowing in general for a nontrivial renormalised ground state energy. 



4 Cosmic string with finite thickness 

We use the metric given by Eq. (|13"D, (PI). Again, the coordinate z will be 
dropped because we work in (2+1) dimensions. The curvature scalar is 

2e 2 

R= T , (31) 

Po 

and we allow for £ ^ |. By means of the ansatz (n = 0, ±1, ±2, . . .) 



i ^m(p)A/ ffi sin(^) , p G [0,p ] 
$ = e mip g~ I TZ = e inLp I V e k po j ' F L ' FUJ (32) 



we arrive at the equations 



d 2 e> 2 - 1/4] e 2 



dp 2 p\ sin 2 (ep/p ) Ap 2 Q 



+ —(l-80 + \ 2 \n in = 0,pG[0,p ] (33) 
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d 2 nV-1/4 
dr 2 r 2 



+ X 2 \TZ 



out 



0, re [r , oo) (34) 



for the radial functions. Both this equations may be solved exactly. Indeed, the 
solution of the radial equation outside of the string (Ell) can be written as 



n 



out 



f n (X)H w (Xr) - m)H+ u (Xr) 



where f n is the Jost function and 

fli(Ar) = ±i\ 



'irXr 



H^ 2 \Xr) 



(35) 



(36) 



are combinations of the Hankel functions. The solution regular at the origin of 
the radial equation (|33D inside the string is 



In . ep ( Xp 
-sin — 

2 Po V e 



n+l/2 



p— n 



cos • 



Po 



with 



a 



and P™ is the Legendre function. These solutions must obey the matching con- 
ditions 

^in(po) = K out {r ) , n' in {p ) = n' out {r ) 

on the surface of the string. From this we get the following formula for the Jost 
function 

n+l 

/ 77 Hill f / A i). . \ 

/n(A) 



_ e i\n(y-\) 



^7^ sin e 



x 



2 ^/cos e \ e J 
^ 1 /(Ar Q )p-«[cos£] + ^)(Ar Q )P^ , [cose] 



e sin e 



Ap c 



Here the primes denote derivatives with respect to the arguments. Taken on the 
imaginary axis this Jost function reads 

n+l 



fn(ik) 



sine 



/ fcpo 

x/cose \ e J 

x |iC(£;ro)P- n [cos e ] + K nv (kr )P- n '[cose 



(37) 



e sin e 



A;p 



Note that this function doesn't have zeros for k e [0, oo), i.e., there are no bound 
states. This can be checked by inspection of eq. (|37j). The Jost function has the 
following asymptotics for large resp. small k: 



fn(ik) 



k— too 



exp{-A;p (— - 1)} 
Po 



fn(ik) 



k^O 



k^O 



Ink . 



(3? 
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Using the formula 

limt n P7 n [cos-] = J n {x) , (39) 
the Minkowskian limit {y — > 1) for the Jost function 

lim/ n (A) = l (40) 

V— *1 

can be checked. 

Proceeding as it was done in the case of the infinitely thin string we obtain 
the following expression for the ground state energy in (2 + 1) dimensions 



E = -M 2 *^X> n [°° dk[k 2 - m 2 ]^'* 

27T n Jm 



n=0 

d_ 

dk 



x-^ln{AT" [f n (ik)H- v (ikR) - f*(ik)H+ v (ikR)]} 



and, using 

after obvious rearrangements we get 



E = -M 2sC -^fd n £dk[k 



m 2ll/2- S 



n=0 

d_ 

dk 

_ M 2s^l g dn r dk[k 2 _ m 2]V 2 -|- In k-"f n (ik)I nv (kR) 

Z7T n= Q J m OK 
271 T^n Jm Ok 



• — In {k~ n [f n {ik)I w (kR) - f n {ik)K nv {kR)] } 



l f n jik) K w {kR) 

f n (ik) I w {kR) 



Now the contribution of the last line is exponentially small for R —>■ oo. The 
contribution of the preceeding line can be written as the sum 

E = El hin + (41) 

of the contribution E^ in ( PD| ) of the infinitely thin string considered in the pre- 
ceeding section and the contribution of the interior structure of the string 

+oo 



£int = _ M 2s^lf dn r dk[k 2 _ m 2 ] l,2-.d hlk n(v-l) fn{ik) (42) 
Z7T n= o J" 171 OK 



After the work done in section 3 it is just this contribution which remains to be 
calculated. 
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5 Approximation of small angle deficit 



To calculate E 1 ^ (f42"D and the corresponding renormalised ground state energy 
the analytic continuation in s to s = must be performed. To this end the 
knowledge of the uniform asymptotic expansion of the Jost function f n (ik) (|37j) 
for n — > oo, k — > oo, n/A; fixed is required. Now, althought this Jost function 
is known explicitely in terms of Bessel functions and the Legendre function, this 
task is not easy to perform. The point is that the asymptotic for a —>■ oo and 
n —>■ oo of the Legendre function P~ n is quite complicated to handle. Therefore 
we restrict ourselfs to the easier tractable case of a small angle deficit resp. mass 
density of the string, i.e., to the case e< 1. Than in the radial equation for the 
interior of the string 



dp 2 



n 2 - 1/4 
~~P^~ 



+ \ 2 -u( P ) \n in = o , pe [o,p ] 



with the potential (using (|31|)) 



U{p) 



n~ 



sin 2 9 



1 + 6 2 2£ 



ep/po) we approximate 



U(p) = U + 0(e 4 



with 



U 



3pg 



[n 2 + 6£ - 1) 



By this, Eq. d43|) can be solved in terms of Bessel functions 



An 



n+1/2 



7T 



-A n p J n (X n p) 



/n 6 (A) 



Uq) and the corresponding Jost function reads for small e 

J n {\ n p Q )H^{\r Q ) 



i—e 2 ( ' 
2 



A_ 

An. 




+ VroPo [KJ' n (KPo)H^(Xr ) - AJ n (A n p )i/,l 1 /(Ar )] } 
On the imaginary axis we get 



' k_ 




In{k n po)K nv (kr Q ) 



ir 



+ J— [k n p o r n (k n p )K nu (kr ) - kp I n (k n p )K' ni/ (kr )} 
Po 



(43) 



(44) 



(45) 
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with r /p = tane/e and k n = ^Jk 2 + fji(n 2 + 6£ — 1). It is easy to verify that 

it obeys the limits fl3§D and (|40|). 

Note that the approximate potential U (fPj) is constant and that Uq may 
take negative values. Therefore bound states could occur. However, due to the 
conical structure of the exteriour part of the space (see v in the index of the 
Hankel functions) and the corresponding relations between the parameters, in 
fact no bound states occur. Correspondingly, it can be shown that the Jost 
function ( f4"5| ) doesnt have zeros for k e [0, oo). 

Now we insert this approximate Jost function into E 1 ^ ([42D and have to 
perform the analytic continuation in s to s = 0. For this reason we use the 
uniform asymptotic expansion of (|45|) which can now be obtained by simply 
inserting the corresponding expansion (^Tj) of the Bessel functions. We define 



f: e ' as m 



-ni>rj(z v )+nri(y) 



24n 



(46) 



with 



rj(z) 



kp 



n 



kr 



;i + , 2 )- 1/2 , 



Vl + z 2 + In 



i + VTT1 2 



, y 



e 2 1 - 6£ 



n- 



In general, the expansion (M) must include all terms up to n~ 2 . But that contri- 
bution doesn't appear (together with n~ 4 , n~ 6 , ...). Next we divide the expression 
for i?Q nt into two parts 

pint j^iint i pint 

simply subtracting and adding In k n ( u ~^ f^' as in fl42|). Here 



(47) 



pint 



—M 



2s cosns 
2tt 



n=0 



dk[k 2 - m 2 } 1/2 - s -^ In k n{v ~ l) f a n s {ik) 



(4£ 



is called the 'asymptotic' part which still requires the analytic continuation to be 
done and 



E 



int 
0,fin 



-1 

2^ 



n=0 



d n dkVk 2 



d_ 

dk 



\nk n ^f n (ik) - kik n ^ f™(ik)} (49) 



is called the 'finite' part. In it the analytic continuation could be performed under 
the sign of the integral and the sum because they are convergent. 

In i?Q n * s (^) the integrand can be expanded in powers of e. Than the inte- 
gration over k and the summation over n can be carried out explicitely using the 
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same method as in section 3. We obtain 



E, 



hit 




7T 

12 



-/r 1 

Q H 24 6 



6£ 



-7T 







12 

ln(l - e~ 2 ^ 



720 



72 



7T 





(50) 



7T 



-(3 2 [3Qi(/3)-3Q 2 (/3) + Q 3 (/9) 



vrQi(/3) 



with the notation /3 = mp and Q a defined by (f2"5"|). We remark that Q a (0) are 
exponentially decreasing for /3 — > oo. 

At this place we can determine the heat kernel coefficients for the thick string. 
We have to consider the asymptotic expansion of for large m. The nonde- 
creasing contributions may be contained only in E^l s ( |4"BD and can be read off 
from Eqs. fl50"D and (|7|). They are 



rjint 



2^0 



mt 



-2<6 2 . 



(51) 



Note that the coefficient Bi/ 2 is zero. This was to be expected because the 
background is smooth enougth not to allow for boundary dependent coeficients 
in this order. Also, the coefficient -B3/2 is zero. This is in the given order in e and 
follows simply from dimensional reasons. In higher orders in e it may be nonzero 
like further coefficients of higher half integer order. 

Now, by means of (41) we have to take the contributions to the heat kernel 
coefficients of the infinitely thin string (^) and that of the 'interior', Eq. (|5l"|) 
together. These coefficients can be compared with that following from the general 
formulas. For instance, from 



P0 K 

B =fdV = 2n f^ S m e -Pdp+- [rdr 

J J € Pq V J 



ttK e 2 7ip 2 



V 



+ 0(e 4 



(52) 



we obtain the boundary dependent contribution (cf. Eq. (p7|) ) and £?Q nt (|51~1) . For 
the coefficient B\ we have 



Bl = i_^, Rd „ 



1 ~6£ 2 f 4 , 
— - — e 7r + U[e . 



(53) 



Now, in the given approximation it holds 



3 V v) 3 
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Therefore, from B\ (|27j ) resp. ( |29| ) of the thin string and B 1 ^ 1 (|5T| ) we get 



7T 



V 



j- 1 27T r, 1 — 

2<e 2 = — + e 2 vr— 
3z/ 3 



6£ 



+ 0(6 4 ). 



(54) 



We note that the Kac term disappeared. The first term in the last line is due to 
the boundary at r = R and the second is the genuine contribution of the thick 
string. It vanishes in case of a conformal coupling. 

To proceed, we define E r Q en by means of (§) using the coefficients B fl52|) and 
B\ (^) in the definition of Eq fl7|). Note that the coefficients with half integer 
numbers are pure boundary dependent contributions resulting from the boundary 
conditions at r = R. 

In fact, the renormalization in E nt is reduced to dropping the nondecreas- 
ing for m — > oo contributions in E 1 ^. Therefore we obtain for the complete 
renormalized ground state energy 



7-ircn 
^0 



1 o2 1 1 ~ 

-8 1 

6^ 24 6 



6£ 



ln(l 



-27r/3\ 



27rpo 

[3Qi(8)-3Q 2 (8) + QM 



- J2 d n 

n=0 



1 1 
+ - 



ttQi(/3) 
-6£ 



720 



6 



(55) 



Here, the notation 



\nf n (ix) - lnf™(ix) 



is introduced. 

Some further work is necesary. In F n the Jost functions have still to be 
expanded for small e. For f n (ik) this can be done using the formula III 



dK p (x) 



n! (x 



dp 



\p=n 



n ~ l 'x\ n KAx) 



E 

1=0 



(n-l)l\ 



for fn{ik) this expansion is a simple task. Finally we obtain 



'(n + If 



x 



n 



x 



x 



n\e=0 



— InKn + I ~ — I n +lK n -i + — I n+ \K, 



nn\ 



n-l 



^ 1=0 
^ 2~ ^ 



i— n 



6 



x(I n+1 Ki + I n Ki +1 ) 



n 



l ~ n InKi 



n-l)l\ 
n 



+ - [K Q n\I n [-) +lnx 



x 



+ 



x- 



1 + -T- 



n 



n>2 



n- 




1 



24n 
16 



i + 4 



l-6£ 



n 



1 + + 



This expression has to be used in i?Q en (|55| 



In writing 



2np 



G(P) + 0(e 4 



(56) 



it is in fact G(j3), a function of one variable, which must be calculated. We did 
this task numerically. After a carefully examination we came to the result 



G(P) = , 

for arbitrary £; the details are given in the appendix. 



(57) 



6 Conclusions 

In the preceeding sections we worked out methods suited for the calculation of the 
ground state energy of a massive scalar field in the background of a cosmic string. 
The main emphasis was on a string of finite thickness. We used the standard 
renormalization scheme, i.e., we calculated and subtracted the contributions of 
the first few heat kernel coefficients. Thereby the normalisation condition, stating 
that the renormalised ground state energy must vanish for a large mass of the 
quantum field, is imposed. 

As a part of these calculations we first considered the infinitely thin string in 
detail. Using explicit formulas we reobtained the known heat kernel coefficients. 
These are the coefficients due to the boundary of a large cylinder at r = R, which 
was introduced in order to render the eigenvalues discrete, and the topological 
Kac term. In the sense of the renormalisation used we got the result that the 
renormalised ground state energy of the pure string, i.e., when removing the outer 
boundary (R — ► oo), is zero. 

Then the same problem is formulated for a string with finite thickness. How- 
ever, the complete calculation suffers still from mathematical difficulties. There- 
fore the approximation of a small mass density e of the string was introduced. 
Having in mind the smallness of e ~ 10~ 3 in cosmological applications, this is at 
once physically motivated. 

We note that in this approximation an alternative calculation should be pos- 
sible, namely the use of a perturbation theory in the mass density as it was done 



for the calculation of the Casimir force between two cosmic strings in . 

In this approximation of a small mass density, in order e 2 , first the heat kernel 
coefficients were calculated. They are checked to coincide with that following 
from general formulas. We remark that the Kac term disappeares and that for 
a conformal coupling (£ = 1/6) there are no counterterms required besides that 
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which follow from the boundary at r = R. From this it is clear that the Kac term 
is due to the singular behaviour of the metric of the thin string at the origin. 
This can be understood from another point of view too. Consider the vacuum 
expectation value of the energy density in the background of the thin string. For 
dimensional reasons it behaves like r~ 2 near the string and, therefore, cannot be 
integrated over r near the origin. Now, if introducing a suitable regularisation, 
zeta functional regularisation for instance, it becomes possible to integrate over 
r. As a result, when removing the regularisation, an additional divergence occurs 
which is just the Kac term. 

After performing the renormalisation we calculated numerically the ground 
state energy in the background of the thick cosmic string in order e 2 of a small 
mass density. The result is zero with the reasonable precision of 10~ 7 indepen- 
dently of the parameter £. Thereby a nontrivial compensation between different 
contributions occured. 

This result that the vacuum of a scalar field is not disturbed by a cosmic string 
is quite remarkable and seems not to have analogues in other configurations. For 
instance, we do not see any symmetry or invariance arguments for this result 
although they should be there. 

Perhaps, there is some relation to the result of Brevik and Jenesen [16j indi- 
cating the absence of particle production in the formation of a cosmic string. 

Further work is necessary. For instance, the result should be extended to the 
(3+1) dimensional case, to higher spin fields and to mass densities which are not 
small. 
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Appendix 



Here we consider the analytical and numerical analysis of the function G{j3) 
defined in Eq. (|5T|) 



G(f>) 



6 H 24 6 



ln(l _ e - 2 ^) - |/5 2 [3QM - 3Q 2 (/3) + Q 3 ((3) 



l-6£ 



nQx(P) - d n dx\]x 2 -f3 



n=0 



;0F n 
dx 



720 



l-6f 
72 



7T / \ 



Obviously we have G^fl)^^ — ► and the domain of interest is the neighborhood 
of zero: (3 ~ 0. For numerical simulations the above formula is more suitable 
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for (3 > 1. In the opposite case, (3 < 1, there exist at first sight a logarithmic 
singularity for small (3 : ln(l — exp(— 2irf3)) ~ ln27r/3. But this singularity is 
canceled with that in the contribution of (n = 0) in the series in fl58|). For 
numerical calculations it is more suitable to cancel this singularity in manifest 
form. For this reason we divide this term into two parts, 



J°° dx^x 2 -(3 2 ^-F = J' dxy/x 2 -/3 2 -^F + J™ dx^x*-[3 2 -^F , 



(59) 



where 



F 
$ 



x 



6 



-IoKo 



x 



x 



X 1 

~3 + 24~r 



The integral 



3 
6x 



dx 



-I K 2 



6£ 



{I K + hK 



(60) 



may be found in manifest form. Thereby we arrive at the following expression 
for the case (3 < 1 



G(/3) 



6 f 24 



l-6£ 
6 



In 



1 - e- 2 ^ 



(61) 



7T 1 — 6£ 

- 5 i3 2 [3Qi(.i3) - 3Q 2 (/3) + Q 3 (/3)] - —^^(/J) 



l-6£ 



with the notation 



Qn(/3) 



P dtt r 



(3 n Jo e 



o2vrt 



(62) 



It is necessary to stress that formula (|BTD is only a different representation 
of G{(3) ( |55D which we made in order to avoid the logarithmic contribution in 
individual terms. Next we observe that the series over n is slowly convergent and 
quite a large number of terms must be taken into account 

Now let us consider some analytical properties of G{(3). This function doesn't 
have a linear term in the expansion for small (3, 



G((3)=G(0) + O({3 2 ) 
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In order to argue this statement let's consider the part in (O) which contains the 
integral 



D(p) = f dx^jx 2 - P 2 -!^+ f°° dx^jx 2 - (3 



dx Ji v dx 
+ E2/ dx^x 2 -(3 2 -^. (63) 

n=l J P ux 

The function $ has the following expansion for small x (using the power series 
expansion of the Bessel functions) 

00 00 

^ k=0 fc=0 

It is only the zeroth term in the first sum, CqIb.x/2, which deliveres a linear 
contribution to the first integral in (|63|) . All other terms contribute higher powers 
in (3. Thus 

f dx^jx 2 -(3 2 ^ = £ dxx^ - ^C Q /3 + 0(f3 2 ) 

with C = (1 - 60/6 - 1/6 = 

In the second integral in (^) we can expand the integrand for f3 2 /x 2 <C 1 
obtain 



dxJx 2 -(3 2< ^- = / dxx^- + 0(13 



dx Ji dx 

In the last term in (|63|) we use the uniform expansion of F n ( pj|) which starts 
from the third power of 1/n. All integrals can be calculated in closed form, their 
expansion starts from the second order in (3. Therefore 

00 roc 1 f)P 00 roo f)fr 

Puting together all three parts we obtain 

D((3) = D(0) + ^ + O((3 2 ) . (64) 
The expansion of the remaining terms in ( |B~1]) contains the same linear term 

const +^/3 + 0(f3 2 ) , 



which is canceled by that in (|64|) . Therefore the expansion of G(f3) starts at least 
from the second power of (3. 



20 



Figure 2: Plots of the integrand N n = 2ny-^F n ((n 2 y 2 + fi 2 ) 1 / 2 ) for a) (3 — and 
b) (3 = 0.4 and n = 1,2,3 

For the numerical analysis of the series in G{j3) first of all we replace the 
integration variable x — > y = \/x 2 — (3 2 /n. Then we obtain 

00 poo J f) 00 roo f) 1 

£ / 2^x 2 -(3 2 —F n (x)dx = / 2m/— F n (Jn*y> + 0>)dy . 
i^xh dx ^Jo dy 

Some first integrands iV n = 2nyJ^F n (y/ri 2 y 2 + /? 2 ) are plotted in Fig.|| for /3 = 
and for (3 = 0.4. As it is seen from the figures, all functions N n (y) have quite 
large variations near the origin and decrease as 1/y 3 for large y. 

For higher n, (n > 3), the Bessel functions entering F n , have been substituted 
by their unifom asymptotic expansions whereby the first 13 terms were taken. 
The error caused by this approximation is smaller than 10 -7 . Than the integral 
and the sums can be carried out explicitely. For this task Maple was used. 

With this, the function G{(3) was calculated for < (3 < 2 and the result 

I G{p) \< 10~ 7 
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was obtained numerically. 
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